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AERONAUTIC SYMBOLS 
1. FUNDAMENTAL AND 'l«a«irt1HB- UNITS 



Time 

Force. 



Symbol 



I 
t 

F 



Metric 



Unit 

i 



meter 

second — 

weight of 1 





English 



Unit 



_ 



foot (or i 
second (or hour)... 
weight of 1 pound. 



_4 * 



horsepower 

miles per hour — 
feet per second.. 



Abbrevia- 
tion 



ft. (or mi.) 
seo. (or hr.) 

lb. • 



lip. 

m.p.h. 
f.p.s. 



Wt Weighty r, Kinematic viscosity 

g, Standard acceleration of gravity = 9.80665 P, Density (mas ,cr unit volarne) 

[• )f „ac 2 andard • • f dry air, 0.12497 kg- at . 

W ° Tj\yM*1 * ■*:>. ... , '.. X5° C. and 700 mm; or 0.0Q2378 lb#t.^ fieo.* 

W8 B '|»; Soecific weight of 'Standard'! airl i.2255 fcg/n# t 



• " •■ "■ 



Elinema 



TO, 



Specific Weight of "standard" airi 1.2^55 *£/ 
I, Moment of inertia = m/t% (Indicate axis of 0.07661 lb./cuifc. Jv-V ; VV- ;" 

■ . ' radius of 'gyratipnfc. by 1 proper subscript.) ">■"-■• % W$UiMKmm 

'• •»• • • » • 



mm 



Coefficient of viscosity 1> ,? , 

/<' • ' '' S. AERODYNAMIC SYMBOLS 



i : ' , 



'Si'- . Area 

Area of 'wing • 



Area of 'wine • v **Jv line) • •!i''-y.?iv i -''« x'S-V: /•• = ' '-■ 

i„ Angle of stabilizer setting (relative to thn it 



i*,' Anglo of setting of wings (relative. ^to -thfuat 



*Jf: 




X, lift, absolute coefficient (7 t = 4 of 10 cm chord, 40 m.p.s'. the corresponding 

^ number is 274,000) 

D, Drag, absolute coefficient C' D = ^ G P , Center-of-pressure coefficient (ratio of distance 

^ of c.p. from leading edge to chord .. fell) 

I>„ Profile drag, absolute coefficient C B ~ -a «, Angle of attack 

*| ti Angle of downwash 

A, Induced drag, absolute coefficient <J D . = «„, Angle of attack, infinite aspect ratio 

V «,, Anglo of attack, induced 

Z>„ Parasite drag, absolute coefficient Cn,-^ a,, Angle of attack, absolute (measured from zero- 

- '-. lift position) 

C, Cross-wind force, absolute coeffic; 

B, irce 



"I.'.' 
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PROPULSION OF A FLAPPING AND OSCILLATING AIRFOIL 

By I. E. Gabkick 



SUMMARY 

Formulas are given for the propelling or drag force 
experienced in a uniform air stream by an airfoil or an 
airfoil-aileron combination, oscillating in any of three 
degrees of freedom: vertical flapping, torsional oscillations 
about a fixed axis parallel to the span, and angular 
oscillations of the aileron about a hinge. 

INTRODUCTION 

It is the object of this paper to investigate theoreti- 
cally the horizontal forces experienced by an airfoil 
or an airfoil-aileron combination in a uniform air 
stream made to execute flapping motion or to perform 
angular oscillations about a fixed axis parallel to the 
span The problem treated is that of an infinite wing, 
or wing and aileron, performing steady sinusoidal oscil- 
lations in any of three degrees of freedom: vertical 
flapping at right angles to the direction of motion, 
oscillations about an arbitrary fixed axis parallel to its 
span, and oscillations of the aileron about a hinge. 

The work of Wagner (reference 1) for calculating 
the distribution of vorticity in the wake of an airfoil 
in nonuniform motion appears as a starting point. A 
vortex wake is generated by the oscillatory motion, 
which in turn affects the, en tire nature of the forces 
experienced by the wing. Beautiful experimental 
checks of Wagner's theory of the manner in which the 
circulation builds up have been obtained by Farren 
and Walker. (Cf. reference 2, ch. 9 for a more detailed 
bibliography.) Birnbaum and Kussner (reference 3) 
have also attacked the problem of obtaining the lift 
forces on an oscillating wing by certain series expan- 
sions that are rather cumbersome to handle. Glauert 
(reference 4) has treated the case of an oscillating air- 
foil and has obtained expressions for the forces and 
moments that check with Wagner. Theodorsen (ref- 
erence 5) has developed compact expressions for the 
lift and moments in the case of an airfoil-aileron com- 
bination of three independent degrees of freedom and 
has applied the results to an analysis of the wing-flutter 
problem. The foregoing references are concerned only 
with the lift forces, not with the horizontal forces; 
however, von Karman and Burgers, who present in 
reference 2 a resume' of the work (to 1934) on non- 



uniform motion, calculate there <the propulsion effect 
on a flapping wing. The present paper makes appli- 
cation of the compact formulas developed by Theo- 
dorsen and of the method outlined by von Karman 
and Burgers to treat the propulsion on a wing oscil- 
lating in three independent degrees of freedom. 

The assumptions underlying the theory are small 
amplitudes in the various degrees of freedom and a 
(infinitely) narrow width of the rectilinear vortex 
wake, as well as the usual assumption of a perfect fluid. 
Quantitative agreement with experimental values, 
which are not very abundant, can hardly be expected 
since the finite width of the wake is important with 
regard to considerations of the resistance; neverthe- 
less the results can be useful for interpreting such ex- 
periments as exist on the so-called "Katzmayr effect" 
(reference 6) and for clearing up certain aerodynamic 
features of the nature of the flight of birds. 1 Experi- 
mental tests on an oscillating and flapping wing are 
being conducted at the present time by the N. A. C. A. 

This paper is not concerned with the problem of 
flutter, which is an instability phenomenon that mani- 
fests itself in certain critical frequency ranges and is 
due to an interaction and feedback of energy because of 
coupling in the various degrees of freedom. (Cf . refer- 
ences 3 and 5.) Profile drag is to be considered as ad- 
ditive to the horizontal forces obtained. 

FORCES AND MOMENTS ON AN OSCILLATING AIRFOIL 

Consider an airfoil represented by the straight line of 
figure I . The airfoil chord is of length 2b and (its mean 
position with b as reference unit length) is assumed to 
extend along the x axis from the leading edge x= — 1 
to the trailing edge x= + 1 . The coordinate x—a repre- 
sents the axis of rotation of the wing, x=c the coordinate 
of the aileron hinge. The airfoil is assumed to undergo 
the following motions with small amplitudes: A vertical 
motion h of the entire wing, positive downward; a 
rotation about x—a of angle of attack a, positive clock- 
wise and measured by the direction of the velocity v 
at infinity and the instantaneous position of the wing; 
an aileron motion about the hinge x=c of angle j3, 

1 It is interesting to observe that the Katzmayr effect occurs in nature also in the 
motion of flsh. See "The Physical Principles of Fish Locomotion," by E. G. Rich- 
ardson, Jour. Exp. Biology, vol. XIII, no. 1, Jan. 1936, pp. 63-74. 
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measured with respect to the undeflected position of 
the wing itself. 




(2) 



Oj Axis of 
rotation. 
I Aileron hinge. 

Trailing edge. 

Figure 1.— Parameters of the airfoil-aileron combination. 

Let us consider sinusoidal oscillations in the various 
degrees of freedom and use the complex-number nota- 
tion 

a- 

^fte'^+V* (l) 

The constants a 0 , /3 0 , and h 0 represent the maximum 
amplitudes in the various degrees of freedom, <p 0 , <p 1} 
and (p 2 are phase angles, and the parameter p determines 
the frequency of the oscillations. By means of the rela- 
tion 

lev 

"J 

an important parameter k is denned, i. e., k=pb/v. It 
will be seen that 2-irjk is the wave length between suc- 
cessive waves in the vortex wake in terms of the half- 
chord b as reference length. 

The following three formulas for the lift and moments 
on an oscillating airfoil of three degrees of freedom are 
due to Theodorsen and are taken from reference 5: 2 



~2w P vbC(k)Q (3) 
= - P b^7r(~ - ajvba + wb 2 (^ +a^a+ T ls v>l3+ T l6 vb$ 
+2T 13 b 2 (i-aM^+2pvb 2 ir(a+^C(k)Q (4) 

M fi =- P b 2 ^T 17 vba+2T n b 2 a+h 2 T ls p-^-vbT l9 0 

-^T s b 2 p-TM~]r-(yvb 2 T l2 C(k)Q (5) 



where 



These equations are to be interpreted as follows: The 
real part of P denotes the lift force (positive downward) 

> The writer wishes to record the fact that in order to establish a check on these 
general relations he has compared them with the widely varying expressions given 
by Wagner, Qlauert, yon Karman and Burgers, and Kiissner in their special cases 
(references 1 to 4). Identical agreement has resulted in all cases, except that in the 
case of Kiissner's formulas a numerical cheek was made sincean analytic check was not 
feasible. The numerical agreement was good except in the case of the wing-afleron 
combination where Kiissner makes some rough approximations. 

A recent paper by Cicala (reference 7) deserves mention. Cicala derives expressions 
for the lift and moment on an oscillating airfoil that seem to agree with the results of 
Theodorsen, although the method is somewhat more involved. The functions de- 
noted by Cicala as X' and X" correspond to 1-F and -Q defined in equation (6). 



associated with the motion given by the real parts of 
(1); i. e., a=a 0 cos (pt+<p 0 ), j8== j3 0 cos (pf+pi), and 
h=ho cos (pi+¥>2). The imaginary part of P denotes 
the lift force associated with the motions a=a 0 sin 
(pt+<po), /3=A, sin (ptf-f-.pi), and h=h 0 sin (pt+w). 
Similarly M a and Mp denote in complex form the mo- 
ments (positive clockwise in fig. • 1) about x=a and 
x=c, respectively, due to the motions (1). (The mean 
value of a or /3 is considered zero. When the mean 
values are different from zero, the forces and moments 
arising from constant values a m and 8 m are to be added.) 
In equations (3), (4), and (5) there occur various 
symbols that have not yet been denned. The T's, 
i. e., Ti, T z , Ti, etc., are constants defined completely 
by the parameters c and a (reference 5, p. 5). For 
reference they are listed in appendix I, where there is 
also given a collection of the symbols employed in the 
notation of this paper. The function GQc) is a useful 
complex function of the parameter k (see (2)) and is 
given by 

C(k)=F(k)+iG(lc) (6) 

wbere 

p _ J 1 (J 1 + Y 0 ) + Y l (Y 1 -J 0 ) 
*- (J 1 + Fo) 2 +(F 1 -J 0 ) 2 



Q-- 



(j l +Y 0 y+(Y t -j 0 r 



Functions J 0 , J u Y 0 , and Y x are standard Bessel func- 
tions of the first and second kinds of argument k. 
Figure 2 and table I, which are taken from reference 
5 (with certain minor changes), illustrate these func- 
tions. 

In what follows we shall be interested only in one 
part of the preceding complex equations. It is an arbi- 
trary matter whether to employ the real or imaginary 
parts. The choice made here is to treat the imaginary 
parts, and we write down for reference the imaginary 
parts of equations (1), (3), (4), and (5): 



a=a 0 sin (pt+<p 0 )\ 
/J=j8 0 sin (pt+tpi) \ 
h—ho sin (ptf+^ 2 )J 



(7) 
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Figt/be 2.— The functions F and —6 against l/k. 
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P=—pb 2 [viraop cos (pt+<p 0 )-- ith^p 2 sin (pt+<p 2 ) + Tbaa 0 p 2 sin (pt+<f>o) 
—vTSoV cos (pt+vJ+TibPop 2 sin (pt+^i)] 

—2irpvb^va 0 sin (pt+<p Q )-\-hoP cos (pi-\-<p2)+b(^—aJa 0 p cos (pf+^o) 
+^-Vo sin (pt+n) + T$ftp cos (pi+p,)] 

— 2xp»6Gj^«a 0 cos (pt+<p 0 ) — hop sin (pf-|-$> 2 )— 6^— a^aoP sin (p<+^ 0 ) 

+^%A> cos (pf+ ft )-|s6flp sin (p«+ft)] (8) 

■M«= — pJ^ttQ— ajvba a p cos (#*+p 0 )— x& 2 Q+a 2 ^p 2 sin (ptf+*>o) 
+ ri 6 « 2 /3 0 sin (ptf+pO + TjePfc/SoP cos (ptf+pi) 
—2T 13 b 2 p 0 p 2 sin (pt +<?!)+ awbhop 2 sin (p<+^ 2 )J 

+2pw6 2 7r^a+^i ? '^»ao sin (pi+*> 0 )+AoP cos (pt+<p 2 ) 
+bQi—a)°QP cos (p«+<Po)+^Vo sin (pt+n) 

+ ^6j3 0 p COS (p« + ^l)J+2p»6 2 7r^a+0(r^»ao COS (pt + <{> 0 ) 

—h 0 p sin a^)aop sin (pt+(p 0 ) + ~vfi 0 cos (p<+<Pi) 

-^Apsin (pt+vi)] (9) 
Mg=— pb^Tnvbaop cos (p«+*> 0 )— 2T n b 2 a {s p 2 sin (p«+Vo) +^« 2 Ti8)3o sin (pi+ft) 

-%vbT n fop cos (p<+^)+^Wsin Cpf+pO + WoP 2 sin (p<+«)] 

— pw6 2 7 , 12 i ? |^»ao sin (pi+^ 0 )+AoP cos (pt+<Pz)+bQ^— a)a 0 p cos (pi+¥> 0 ) 

+ ^H>/3„ sin (p<+¥>i)+^>A,P cos (p<+^)J— p^Tja^ao cos (pt+<p 0 ) 

—hp sin (pif+*5 2 )— &Q— ajjaop sin (p«+^ 0 )+^/3o cos (p<+¥>0 

-^bfop sin (10) 

In addition to these equations we will need the expression for the force on the aileron. This equation is obtained 
in complex form as (use formulas on pp. 5-8, reference 5) 

P,=- P b 2 (^-vT i a~Tji+bT 9 a-~vT s $-^T 2 ^ 

-2pvbT w C{k)Q 
And the imaginary part is 

Pp=— pb 2 ^— vTiOoP cos (pf+<p 0 ) + 2V&oP 2 sin ipi+<p 2 ) — bT^p 2 sin (pt+<p 0 ) 

-^vT s Pop cos (pt+ Vl ) +^T 2 p 0 p 2 sin (pf+^o]-2p&«Vl^ 2 [|(l-e)a<oP cos (pt+<p 0 ) 

+~Vt/1—c?Po sin (pt+pi) +^.(1— c)T w p 0 p cos (pi+*>i)J— 2pvbT m F^va 0 sin (pt+<p 0 ) 

+h 0 p cos (pt+<p 2 )+bQ^— a^aop cos (pt+<po)+-^vfa sin (pf+pi) 

+^&/?oP cos (p*+pi)J— 2pw6Z , 2oG^»«o cos (pt+<p Q )— h 0 p sin (pt+<p 2 ) 

~-bQ,-a)a 0 p sin (pJ+W+^A. cos (p«+v>i) — -^ftp sin (p<+*>i)J (11) 
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Following the method of von Karman and Burgers 
(reference 2), the average horizontal force will be 
determined in two ways: (1) by the energy formula 
given in equation (12), and (2) by the force formula 
given in equation (13). The agreement of the results 
of the two methods will thus furnish a check on the 
work. 

ENEKGY FORMULA 

W^E+P x v (12) 
where W represents the average work done in unit 
time in maintaining the oscillations (7) against the 
forces and moments (8), (9), and (10) ; E represents the 
average increase in kinetic energy in unit time in the 
vortex wake and; P x v denotes average work done in 
unit time by the propulsive force P x ? 

FOHCE FORMULA 

P x =7rpS 2 +aP+pPi, (13) 

where P x is the propelling force; a and 0 are given in 
(7); P in (8), and Pg in (11); S is obtained from the 

relation 8=)im~ y-Jx+i where y is the vorticity 

distribution. The value of S is finite, since y is infinite 

in the order of at the leading edge x= — 1, and is 

given in equation (25) and derived_in appendix II.* 

We proceed first to evaluate W in equation (12). 
The instantaneous rate at which work is done in main- 
taining the oscillations is , 

W= - (Ph+M a a+]\dp$) 

For the average work done in unit time we have 

2tt 

W=~^jjph+M a a+M p $)dt (14) 

On employing equations (7) to (10) and performing the 
indicated integrations, we obtain after some lengthy 
but elementary reductions 

W=Trpb 2 ^(B 1 h 0 2 +B 2 a 0 2 +B 3 (3 0 2 +2B i a 0 h 0 +2B 5 (3A 



where 

Bi=F 



+25 6 a 0 /3 0 ) 



(15) 



B 2 

£ 3 



="{KH~KX<H+f] 

i 2 | 1m i Tu/liij-, . T i0 G\"~\ 
= 6 L~4^ 2+ ^rV2^ F+ Vi;J 



3 When the energy released in the wake in unit tim e is less than the work required 
in unit time to maintain the oscillations, i. e., EKW, then Px is positive and is a true 
propelling force. When E> W, then P* is negative and denotes not propulsion but 
resistance or drag. 

* Formula (13) is obtained by a slight extension of the method of reference 2, pp. 
305-306. The "suction" force jtpS 2 arising from the infinite vorticity at the Ieadingedge 
is explained in reference 2 (pp. 52 and 306) along lines laid down by Grammel and 
Cisotti. (See also reference 8, pp. 135 and 203.) The fact that this taflnity occurs 
implies that the ideal flow for an infinitely thin wing is unrealizable. We are regard- 
ing this ease, however, as a limiting one of a wing that is rounded and smooth at the 
leading edge and sharp at the trailing edge. 



-(f-G^sin (<p 2 -vo)] 



B *~ 2lL4^~V2- a ;^r+V2^ _ ~2^~ a ) F 

+( ? ^-v a >] sinfe -^ ) } 

In order to calculate Sin equation (12), we need the 
expression for the vorticity in the wake. The magni- 
tude of the vorticity in the wake is given in complex 
form by 

U=U 0 e^e v • ) ( 16) 
where Uge 1 * is a complex quantity determined in (19). 
(Cf. reference 5, p. 8, in which x instead of x—l is used 
in the exponent.) From the definition of the circula- 
tion about the airfoil as the integral of the vorticity in 
the wake we have in complex form, 

T^J'Udx^^UaeW (17) 

Also from reference 5, equation (8), the condition for 
smooth flow at the trailing edge leads to the relation 

+ b^^Q (18) 
Combining (17) and (18) we may write 

r=2^-p( A /|±|-i)^ 

On equating coefficients of e ipt on both sides of this 
relation and solving for the quantity TJtffr (for the 
evaluation of the definite integral in terms of Bessel 
functions see reference 5, p. 8), we obtain 

U^= -4(4+iB) (J+iK)e~ ik 



where, 



(19) 



and 



B 



D= (J, + Y a f+ (Fi-Jo) 2 , J 2 +K 2 =± 

va 0 cos <p 0 —h 0 p sin <p 2 —b(^~—aj(x 0 p sin <p 0 
+— vfio cos n—Zpbbp sin 

=va 0 sin tpo+h 0 p cos <p 2 -\-b(^~—aJa 0 p cos <p a 

T T 
+— °»i8o sin <Pi+~bfi g p cos Vt 
7T 6ir 



(20) 
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When the imaginary part of U is denoted by y, which 
is the only part of interest, the vorticity in the wake is 
given by 

y=Aa cos kx+B 0 sin kx (21) 

where 

1 

-r 



U 0 = {BK-AJ) sin pt- (AK+BJ) cos pt 



±B 0 =(BK-AJ) cos pt+ (AK+BJ) sin pt 

The induced vertical velocity at a great distance x 
downstream is 



\ dx' '=t:(Ac, sin kx—B 0 cos kx) 
x—x 2 V u / 



Tbe difference in potential at points of the x axis in the 
wake is 

<!>2—4>i=bJ'ydx=j : (A 0 sin kx—B 0 cos kx) 

and the kinetic energy in the wake (per unit length) at 
a point x along the surface of discontinuity far from the 
airfoil is 6 



E^PMnfa-h) (22) 
=^p^(^l 0 sin kx—B 0 cos kx) 2 



=4f^[(BK-AJ) cos (pt+kc) 

+ (AK+BJ) sin (pt+kx)f 

The mean value of E x with respect to time is inde- 
pendent of x and is given by 



And, finally, the average value of the increase in energy 
in the field in unit time is 



or also 



E=vE 1 J^(A'+B") 



E=Trpb^[C l h s ?+CW+W+2C i aJi 0 
+2C s AA>+2C 6 «oA>] 



(23) 



where 



■kD 



G_ Js[p + G~ a )] 

«=s[@*+(ay] 

^ 4= ^i^[~l sin fe— <Po)+(|— a^cos fe-Wj 

Equation (12) now defines P^a and hence P x . We 
have 

or 

P x =ir P bp 2 [A 1 h 0 2 +A 2 a Q 2 +A 3 (3 () 2 +2A i a 0 h l) 

+2A 5 0 o h o +2A 6 a o p o ] (24) 
where from equations (15) and (23) 

A X =B X -C X , A 2 =B 2 ~0 2 , etc. 



We shall now proceed to the direct calculation of P x from (13). The value of S is derived in appendix II 
and in complex form is given by 

S= ^2C(k)Q-ba~Jl=?vfi+T 4 b$] 
Again we shall use only the imaginary part, of this expression which is 



•A 



S= (M sin pi+N cos pt) 



(25) 



where 



M=2/'jjto { , cos <p 0 -h 0 p sin ^ b^—a^aop sin <Po+~f/3„ cos <p l — ^bL3 0 p sin <p,J 
— 2G^»«o sin <Po+lhP cos <p 2 + &^— a^oop cos (po+^'^/Jo sin ^+^6^ cos 



+6a 0 ^ sin -yl— c 2 »/3 0 cos <p { — -bfop sin ^, 



fi The expression ^/>w*(02—0i) is actually equal to ^p^~ taken along the 'surface of discontinuity (tbe x axis) where <t>=<t>2—<lti and^=w». The latter expression is 



equal to -^p |~ (gj) + (jy~z) ] taken over a proper space interval, i. e., represents the kinetic energy in a certain volume. 
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and 

A/=2i^»a 0 sin vo+hop cos (p2+b(^-~a^a 0 p cos <po+^vf} 0 sin pi+^bfap cos 

+26 r ^a 0 cos <p 0 —h 0 p sin <p 2 —b(^—a^ceoP sin <p 0 -\-—v^ cos <pi—-'^b$<>P sin pij 

—baop cos <p 0 — Vl — c 2 »j3 0 sin — 4 bp 0 Va cos ^ 
The mean value of wpS 2 with respect to time is 

This expression becomes, after a considerable number of terms cancel, 

irpS 1 2 =Trpbp 2 (a 1 ho 2 +a 2 a o 2 +a 3 ^ o 2 +2a i aJio+2a B pJio+2a 6 a < )0 a ) 



where 



ai =F 2 +G 2 



a 2 =b 2 - 



(^+4^+(i-a) 2 ] + |-(i-a)F-i( ? j 



a 5 =6 



(i* 12 +£ 2 )£— | sin fe;— <po)+Q"-o^) cos fe— <po)~j— | cos (<p 2 — ^ 0 ) + | sin (v 2 — <Po) 

(F 2 +67 2 )[-?|sin te-^+l^cos fa- yi )]+|[ 2 ^ / 7 <;2 sin (%-^)+^cos 

+f[~]p^ cos (<p 2 — (fi)-^sin (^ 2 ~^i)J| 

a 6 =6 2 ((^+^)(g-l+^-)] cos ^-^+[v(l- a )-Sl sin <«-*>) 

, 2*T-2yF=? , 2V1 \ , x pyi^cYl \ . Ti+Tull . , , 

+ 2L— jF"+^(2- a j+2^J Sm ^-W+L^,— (2-«j+-V L Jl C0S 



+ 2-k° I sin ^o)-^ cos (W'-**) 

We proceed to calculate the average values of the 
terms aP and pPp in (13) by employing equations (7), 
(8), and (11). There results 

^P=Trpbp 2 (b 2 a 0 !! +2b l aA+2b i a^ 0 ) (27) 

where 



'2 ¥ 



&4==||^0+^ cos 0p 2 — ^O+^sin (<p 2 — *?o)j| 
, b\( Ti FT xa .O T n \ , s 



Also 



where 



i8P0=7rp62> 2 (c 3 /3 o 2 +2c 6 Mo+2c 6 a o 0o) 



-c 2 

A 2 7T 2 



g Tii^ol 
' 2tt 2 J 



(26) 



(28) 



C3 ~°L 2^ 2 A 2 



£20 

7r 



|Jisin (<pi — spo)} 
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Finally from (13) the average propulsive force is 



P x =irpbp 2 [a l h 0 2 + (a 2 +&2)ao 2 + (a 3 +c 3 ) A, 2 
+ 2 (a 4 + a 0 A 0 + 2 (a 5 + c 5 ) ftA) 
+2(a 6 +6 6 +c 6 )ao/3 0 ] 



(29) 



(30) 



In order that equations (24) and (29) agree we must 
have that 

A 2 =a2~\-b2 
A 3 =a 3 -j-c 3 

A 5 =^a s +c h 

A s =a a +b 0 +c 6 

Each of these relations may be reduced to an identity, 
e. g., consider A t and <h. From (15), (23), and (26) 

In order that Ai = a x the following relation must hold 

M +*+a» (31) 

To show that this is true note that 

F _ J 1 (J 1 + F 0 ) + r 1 (F 1 -Jo) J^+F^+^Fq-JoF. 
(e/i + F 0 ) 3 +(F 1 - t / 0 ) 2 D 



F 2 +G 2 =(F+iG) (F-iG) = 



J 1 2 +Y 1 2 



D 

(cf. reference 5, p. 8) and from a well-known property 
of the Bessel functions, 

2 



e/lFo — e/oFi = 



irk 



Hence equation (31) follows. 

By the use of the relation (31) and the definitions of 
the various T's given in the appendix, it can be veri- 
fied that the remaining relations in (30) are also iden- 
tities. 

It may be of interest to consider the special cases of 
one degree of freedom. Let the motion of the wing 
consist only of the vertical motion h at right angles to 
the direction of flight, i. e., flapping motion. The pro- 
pelling force is then 6 



P x =TC P bp 2 h 0 2 {F 2 +G 2 ) 

Consider in this case the ratio 

P x v energy of propulsion 



(32) 



W 



total energy 
F 2 +G 2 



F 



(33) 



This function, shown in figure 3, represents the theo- 
retical efficiency of the flapping motion (unity = 100 per- 
il This result agrees with the formula of von Karman and Burgers (reference 2, p. 
306) . The expressions of reference 2 denoted by 

h=l+Ai-\ Ai (Q-S)+ X As (P-C) 

&i=Aj-X A, (Q-S)-X A, (P-C) 
reduce in our notation simply to IF and 20, respectively. 



cent). It is observed that a propelling force exists in 
the entire range of l[k, the efficiency being 50 percent 
for irifinitely rapid oscillations and 100 percent for 
infinitely slow flapping. 



1.00 




Figure 3.— The ratio of energy of propulsion to the energy required to maintain the 
oscillations ^pT^ 8 function of Vk for the ease of pure flapping. 

For the special case of angular oscillations about a 
alone (A=0, /3=0) the horizontal force is 

p.^ipwjy^+ffjrp+Q-o)*] 

+KH-p-(Hfl < M > 




Figure 4— The ratio of the energy dissipated in the wake to the energy required to 
maintain the oscillations (SIW) as a function of Uk for the case of pure angular 
oscillations about x=a. 

In figure 4 there is shown the ratio E/W, in this case 
for several positions of the axis of rotation. These 
curves give the ratio of the energy per unit time released 
in the wake to the work per unit time required to 
maintain the oscillations. In the range of values 
0<.E/W<1, Px is positive and denotes a thrust or 
propelling force; for other values it is negative and 
denotes a drag force. 



Langley Memorial Aeronautical Laboratory, 
National Advisory Committee ioe Aeronautics, 
Langley Field, Va., May 4, 1936. 



APPENDIX I 



NOTATION 

a, angle of attack (fig. 1). 
/?, aileron angle (fig. 1). 

h, vertical distance (fig. 1). 

da .. d 2 a • , 
a=z dt' a== d¥' etc " 

«o> ft) ^0) Poi <Pi, <P2, amplitudes and phase angles 
of the oscillations (equation (7)). 

b, half chord, used as reference length. 

x, coordinate in direction of airfoil chord. 
t, time. 

v, velocity of the general motion in direction of 
x axis. 

p, 2x times the frequency of the oscillations. 

k, reduced frequency (equation (2)). The wave 

length between successive waves in the vortex 

wake is 2irbjk. 
a, coordinate of axis of rotation (fig. 1). 

c, coordinate of aileron hinge (fig. 1). 

i, imaginary unit y~ 1. 

e, base of natural logarithms. 
p, mass density of air. 

P, lift force on airfoil (+downward in fig. 1). 
M a , moment on airfoil about a (+ clockwise in fig. 1). 
M,j, moment on aileron about c (+ clockwise in 
fig. 1). 

Pp, lift force on aileron (+ downward in fig. 1). 
OQc), F, 6, J 0 , Ji, T 0 , Ti, J, K, D, Bessel functions of 
the argument k. (Cf. equations (6) and 
_ (19), fig. 2, and table I.) 
W, average work done in unit time in maintaining 
_ the oscillations. 

E, average increase of kinetic energy in the wake 
_ in unit time. 

P x , average force in the direction of the x axis 
(+propulsion, —drag). 
A l ...A 6 ,B 1 ...B fi ,Ci...Gfi, a u ... a 6 , b 2 , b i; b B , c 3 , c 8 , c 6 ,coefficients. 

Q, defined by equation (18). 
A, B, defined by equation (20). 
Ao, Bo, defined by equation (21). 
M,-N, defined by equation (25). 

17, distribution of vorticity in the wake in complex 
form (equation (16)). 
8 



Z7 0 e i *', coefncient'of U, given in equation (19). 
y, imaginary part of U. 

S, defined by equation (13); see also appendix II. 
r, circulation about the airfoil, defined by equa- 
tion (17). 

DEFINITIONS OF THE T's 

Tt=:-^(2+c 2 )^l^?+c cos- 1 c 

T 2 =c(1-c 2 )-(1+c 2 )t/T^? cos- 1 c+c(cos- ! c) 2 
[r 2 =r 4 (T n +T 12 )] 

r 3 =-^(l-c 2 )(5c 2 +4)+^c(7+2c 2 )VT ::r 3 2 cos" 1 c 

-(i+c^cos- 1 c) 2 
Ti—c^l— c 2 — cos" 1 c 

T B = — (l-c?)+2c- } pL = ? cos" 1 c-tcos- 1 c) 2 
T 6 =T 2 

7 1 7 =^c(7+2c 2 )VT ir ?-(|+c 2 ) cos- 1 e 

T s = +2c 2 )tJi^c 2 +c cos" 1 c= -|(1 -c 2 )*-cy 4 

T9=|j(l-c 2 )s +aT 4 ] 

T m = yi^ 2 +CQS- 1 c 
2 7 11 =(2-c)yi : ^c 2 +(l-2c) cos- 1 e 
T 12 =(2+c)i/l^?-~(l+2c) cos- 1 e 
[T 12 -T„=2T 4 ] 

14 — 16 2 fflC 
2\,=2 , 1 -2' 8 -(c-a)2 , 4+|2',i 

[r 16 +r ir = -(|-«)r4+|r„] 

T 17 =-2T 9 -r 1 +(a-|)T4 

Txs=Ts — T±Tio 
Ti§=TiTii 

T 20 = — yT^+cos- 1 c 

mo=T m -2^T=?} 



APPENDIX 11 



EVALUATION OF S (EQUATION (25)) 

From reference 5 (p. 7) we nave that the condition 
for smooth flow at the trailing edge is obtained from the 
equation 

where the <p's are as follows (a ± sign is to be prefixed 
to each <p, + for the upper surface, — for the lower 
surface) : 

1 C" 1 T1 , 

dx 2*Ji V l-x 2 x a -x Uax ° 

(p a =vab-Jl — x 2 

(Ph=hb^l—x 2 

ip- a = «& 2 Q X — a) V 1 — x 2 



Vfi—-vpb[^l—x 2 cos _1 c— (a;— c) log iV] 



1 



^-=2^i86 2 [ Vl -c 2 Vl -x 2 +(x-2c) Vl -x 2 cqb^c 
-(x-c) 2 log iV] 



where 



2V< 



1— CX— -y/l — C 2 -y/l — X 2 



Condition (1) leads to the relation (cf. (18)) 

The leading-edge vorticity may be written as 



(2) 



2S 



On substituting for the p's, making use of relation (2) 
and of equation XI, reference 5, which is 



G{k) 



i Vxp 2 — 1 



X 



Xp + 1 

vV-i 



e~ ikXa dx 0 



there results 

S=^[2C(k)Q-ba- vp+^bp] 
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TABLE I. — VALUES OF THE BESSEL FUNCTIONS 
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k 


1/k 


Ja 


Jl 




Yi 


D 


F 


—Q 


CO 


0 










0 


0. 5000 


0 


10 


Ho 


— 0. 2459 


0.0435 


0. 0557 


6. 2490 


.2548 


.5006 


.0124 


6 


H 


.1506 


— . 2767 


— . 2882 


— . 1750 


. 4251 


.5017 


.0206 


4 


Yi 


— . 3971 


— . 0660 


—.0169 


.3979 


.6389 


.5037 


.0305 


2 


H 


.2239 


.6767 


.5104 


— . 1070 


1. 2913 


.5129 


.0577 


I 


l 


.7652 


.4401 


.0883 


.7812 


2. 6706 


.5394 


.1003 


.8 


W 


.8463 


. 3088 


-. 0888 


.9780 


3. 4076 


.5541 


.1166 


.6 


m 


. 9120 


.2867 


— . 3085 


— L'2604 


4.7198 


.5788 


. 1378 




2 


. 9385 


.2423 


— . 4446 


—1.4714 


0. 8486 


.6979 


.1507 


■4 


m 


.9604 


. I960 


— . 6060 


—1.7808 


7.6823 


.6250 


. 1650 


.3 


m 


.9776 


.1483 


— . 8072 


—2.2929 


11.130 


.0650 


. 1793 


.2 


5 


. 9900 


.0995 


—1.0810 


—3*3235 


19. 570 


.7276 


. 1886 


. 1 


10 


. 9975 


.0499 


—1.5342 


- 6. 460 


57. 810 


. 8320 


. 1723 


.05 


20 


.9994 


. 0250 


— 1.979 


-12.8 


194. 26 


.9090 


. 1305 


.028 


40 


.9999 


.0125 


—2. 430 


—25. 0 


713.4 


.9545 


.0872 


.01 


100 


1.000 


.0050 


—3.006 


—63.7 


4195 


.9824 


.0482 


0 


co 


1.000 


0 


CO 


. — CO 


CO 


1.000 


0 
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Positive directions of axes and angles (forces and moments) are shown bv arrowt 



Absolute coefficients of moment 
(rolling) (pitching) 



a- 



qcS 




'-■ ■ Axis : 


Force 
(parallel 
to axis) 
symbol 


Moment about axis 

... : ': ; 


Angle 


Velocities 


Designation 


Sym- 
bol 


Designation 


Sym- 
bol 


Positive 
direction >; 

; #if|§l|p 


Designa- 
tion 


Sym- 
bol 


Linear 
(compo- 
nent along 
axis) 


Angular 


Longitudinal... 
Lateral.. 


"Jt" 
Y 
Z 


X 

Y , 
Z . 


Rolling .. 

Pitching.... 
Yawing . 


L 
M 
N 


Y >Z 

z — >x 

X >Y 


Roll 

Pitch 

Yaw 


0 


u 
w 


P 
1 
r 


Normal 





D, Diameter 

p, Geometric pitch 

p/D, Pitch ratio 

V', Inflow velocity 

V„ Slipstream velocity 

T, 

Q, 



. control surface (relative to neutral 
C n = ~ position), 5. (Indicate surface by proper subscript.) 

(yawing) 

4. PROPELLER SYMBOLS 

P, Power, absolute coefficient C F =— 
€„, Speed-power coefficient = ^ 



Thrust, absolute coefficient C r ~---r«-, 

pn 2 D* 

Torque, absolute coefficient C Q = -4L« 



leYl 

v, Efficiency 
n, Revolutions per second, r.p.s. 

Effective helix angle = tan -1 { 

\27rr»/ 



5. NUMERICAL RELATIONS 



1 bp. = 76.04 kg-m/s = 550 ft-lb./sec. 
1 metric horsepower = 1 .0132 hp. 
1 m.p.h. -0.4470 m.p.s. 
I m.p.s. = 2.2369 m.p.h. 



1 lb. = 0.4536 kg. . 

1 kg = 2.2046 lb. 

1 mi. = 1,609.35 m = 5,280 ft. 

lm=3.2808ft. 



